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Using a path integral approach, we derive an analytical solution of a nonlinear 
and singular Langevin equation, which has been introduced previously by P.-G. de 
Gennes as a simple phenomenological model for the stick-slip motion of a solid object 
on a vibrating horizontal surface. We show that the optimal (or most probable) paths 
of this model can be divided into two classes of paths, which correspond physically 
to a sliding or slip motion, where the object moves with a non-zero velocity over 
the underlying surface, and a stick-slip motion, where the object is stuck to the 
surface for a finite time. These two kinds of basic motions underlie the behavior of 
many more complicated systems with solid/solid friction and appear naturally in de 
Gennes' model in the path integral framework. 



I. INTRODUCTION 

We study an old but still only very partially understood problem: the dynamics of a solid 
object moving over a solid surface. In practice this is a very complicated and as yet unsolved 
problem, although there is a wealth of experiments, since the general problem is very old 
and ubiquitous in nature, ranging from geology to physics and biology. The basic difficulty 
lies in the very complex nature and behavior of the solid/solid interface, which leads to a 
complicated stick-slip motion of the object [lj]. 

Following P.-G. de Gennes we study in detail one of the simplest phenomenological mod- 
els, far from those of most practical interest, but as a starting point to develop a new 
theoretical approach to describe basic aspects of the above mentioned problem. Ignoring all 
details of the solid/solid interfacial layer, de Gennes proposed a simple Langevin equation 
for the velocity v(t) of a solid object of mass m on a horizontal vibrating surface 0, 3]: 



mv(t) + av{t) + a[v(t))A F = f (t). 



(1) 
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In this equation two kinds of friction between the object and the surface, over which it 
moves, appear: a) a dynamic friction (sometimes called kinetic friction), which is taken 
proportional to v as in the Stokes friction term in fluids and characterized by the dynamical 
friction coefficient a; b) a static friction (sometimes called dry friction), which is given by 
the (t(v)Af term. Here, the function a(v) is the sign function of the object's velocity v, i.e., 
a(v) = +1,0, —1 for v > 0, = 0, < 0, respectively, and A^ is the coefficient (strength) of the 
static friction. The a(v)Ap term represents a nonlinearity and, in fact, a singularity in the 
Langevin equation (pQ), since a(v) is discontinuous at v — 0. Physically, this term ensures 
that the solid object is subject to a static friction, which is equal to A^ and acts always, 
via er(f), opposite to the direction of motion of the object. Both friction coefficients a and 
Ap are here assumed to be constant, which implies that the object moves over an isotropic 
surface. 

In addition, the motion of the object is driven by externally applied one-dimensional 
vibrations of the underlying surface, represented by an external noise £(£), which has the 
properties of Gaussian white noise: 



with noise strength D. The complicated solid/solid interface is therefore replaced by a static 
and a dynamic friction term, and randomness is externally induced by Gaussian noise. 

In this article we use a path integral approach to study the properties of the nonlinear 
de Gennes' model Eq. ([1]). While de Gennes has used a Fokker-Planck approach to obtain 
approximate results for the transition probability (defined in the next section), the path 
integral approach provides a more dynamical picture of the statistical properties of the 
object, on the basis of the most probable or optimal paths in the velocity-time plane. Using 
these optimal paths, we obtain an analytical solution for the transition probability in the 
saddle-point approximation for small values of D. As one of our main results we show that 
the optimal paths of Eq. (CQ) can be divided into two classes of paths, which correspond 
physically to slip motion, where the object moves with a velocity v ^ over the underlying 
surface, and stick-slip motion, where the object is stuck to the surface with v = for a 
finite time. The existence of these two kinds of basic motions is a basic element of almost 
all dynamical systems with solid/solid friction and appear naturally in de Gennes' model in 




(2) 
(3) 
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the path integral framework. 

In the following we present a detailed account of the path-integral approach to nonlinear 
stochastic systems. We analyze the structure of the optimal paths of Eq. (CQ), and derive an 
analytical expression for the transition probability, defined in the next section. 



II. PATH INTEGRAL APPROACH 



The transition probability or propagator f(v,t\vo,to) gives the probability to find the 
object with velocity v at time t, given that it had a velocity vo at the initial time to. Using 
very many paths generated by the Gaussian white noise the transition probability can 
be obtained empirically from many realizations of the de Gennes equation ([1]) for fixed initial 
and final conditions. In the asymptotic time limit t — ► oo the transition probability converges 
to a stationary distribution p(v), which can be derived from Eq. (OQ). For, introducing an 
effective potential 



v 2 



U(v) = — + MA, (4) 

with the characteristic inertial time r m = m/a and A = Ap/m, Eq. ([I]) takes the form of 
Brownian motion in the nonlinear potential U(v ): 

v(t) = -U'{v)+i(t)/m. (5) 

Due to the confining character of U(v ) a stationary distribution of the velocity coordinate 
exists and can be calculated from Eq. §5§ using standard methods 4|. The result is 

p(v) = Ne^ u{v \ (6) 

where 7 = m 2 /D and iV is a normalization constant. Clearly, the stationary distribution 
p(v) is symmetric under a change of sign of v. In fact, also the propagator f(v, t\vo, t ) 
has to be symmetric under the change Vq — > —Vq and v — > —v. This forward/backward 
symmetry is physically due to the fact that the surface is assumed to be isotropic and the 
applied noise is symmetric, so that no bias in forward or backward direction is induced. As 
a consequence, all the statistical properties of the velocity also have this forward/backward 
symmetry. 

In order to obtain the transition probability we use a path integral approach, which was 
introduced into Statistical Mechanics by Onsager and Machlup [s|,[Jto study the fluctuations 
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of a system in thermal equilibrium. It was generalized to fluctuations in systems in a NESS 



in Refs. 



SB 
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ll| . In this approach the transition probability f(v,t\vo,t ) is formally 



expressed as a path integral, i.e., as an integral over all paths leading from the initial state 
(t>o, ^o) to the final state (v, t). For the dynamics of Eq. ([[]), this path integral is given by 



f(v,t\v ,to) = ! {V,t) J[v]e^ A ^Vv, 
J(yo,t ) 



(7) 



where A[v, v] is a functional of v (s), 

A[v,v] 



to 



£(v(s),v(s))ds, 



(8) 



which is usually referred to as the action associated with the path v(s). Here, £ is the 
Lagrangian 



C{v,v) = \{v + U\v)f. 



(9) 



In Eq. ([7]) the integral / Dv denotes an integral over all paths v(s) from (t>o,£o) to (v,t). 
The Jacobian J[v] originates from the transformation — > v(t) and is a functional of v(s) 
due to the nonlinearity of the force —U'(v) in Eq. (jSJ) [12|, Il3| : 



n i k r u"(v(s))ds 

J[v\ = e 2Jt o 



(10) 



We evaluate the path integral Eq. (j7|) in the saddle-point approximation, which proceeds 
as follows (cf. }14|). For large 7 the dominant contribution to the path integral is due 
to a particular path v*(s) that maximizes the exponent in Eq. (JTj), or equivalently, which 
minimizes the action v4[?),t>]: 



5A[v*,v* 



0. 



(11) 



This condition yields an Euler-Lagrange (EL) equation 

d dC dC 

0, 



(12) 



dt dv* dv* 

for the path v*(s), which is the path with highest probability among all paths connecting 
(t>o,to) and (v,t), i.e., it is the most probable or optimal path. We can then expand the 
action A [?),?;] in the neighborhood of the optimal path using 



v(s) = v*(s) + z(s), 



(13) 
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where z(t) is the deviation from the optimal path. Clearly, the boundary conditions for z(s) 
are 2 (to) = z(t) = 0. Expanding the action around v*(s) yields 



A[v,v] = A[v*,v*] + I ds 



* . 5 A 

Sv(s) 



v* 



1 r l r l 
+- ds ds' 



5 2 A 



z(s)z(s') + .. 



(14) 



2 Jto Jto 5v(s)5v(s') 
Here, the linear term vanishes due to Eq. (ITTj) and, using Eqs. (JSJ) and the second order 



term can be calculated as 



52 A 



2 Jt Jt 5v(s)8v(s') 



z(s)z(s 



ds 



to 



z(s 



fl(v*(s))z(s) 



(15) 



with 



n{v) = u"(vy + u'(v)u"\v) 



The leading orders in the expansion of the action are thus 



(16) 



A[v,v] = A[v*,v*]+ / ds z(s) 2 +Q(v*(s))z{s) 



to 



(17) 



Substituting only the zeroth order term of this expansion into the path integral Eq. ([Tj) then 
yields the saddle-point approximation of the transition probability f(v, t\vo, t ): 



f(v,t\v ,t ) = J[v*}e-^^, 



(18) 



valid for large 7. Keeping, in addition, the second order term in Eq. (fTTl) yields the corrected 
form 



/MKto) = J[v*)e-^*>^F[v% 



where the fluctuation factor F[v*} is determined by the path integral 



(19) 



F[v<) 



-7 f* ds\z(s) 2 +n(v*(s))z(s) 2 } 

e *o L J 



Vz. 



(0,t ) 



(20) 



The analytic calculation of the fluctuation factor for the nonlinear potential Eq. (j3J) is beyond 
the scope of this article. In the following we focus on the properties of the optimal paths 
and neglect the second order term in the expansion Eq. (jl7jl . 
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III. SOLUTION OF THE EULER-LAGRANGE EQUATION 

For the Lagrangian Eq. (Q the EL-equation assumes the explicit form 

v* A 

v* g — a ( v *) — = 0. (21) 

A complete picture of the properties of the optimal paths and, on the basis of Eq. (fTBj) . 
of the transition probability f(v,t\vo,t Q ) is obtained by solving Eq. (!2~TT) under the given 
boundary conditions, which are fixed initial and final velocities (vq, to) and (v , t), respectively. 
Eq (1211) consists in fact of two different equations, namely one for positive v*, in which case 
a(v*) = +1 and one for negative v*, where cr(t>*) = —1. Each of these two equations are 
straightforward to solve, if the boundary conditions are such that v* is always positive or 
negative, i.e., if v*(s) remains entirely on either the upper (v > 0) or the lower half (v < 0) 
of the velocity-time (v-s) plane. In that case one finds the solutions (dropping the * for 
optimal path in the following) 

v±(s) = B±e s/Tm + C±e- s/Tm T Ar m , (22) 

where the constants B± and C± are determined by the boundary conditions (cf. Appendix lAl) 
and + refers to the upper half plane and — to the lower, respectively. 

The forward/backward symmetry of the velocity statistics (discussed below Eq. (jB])) im- 
plies that the basic solutions v + (s) and u_(s) are symmetric with respect to the v — axis: 
the path v + (s) between (vo,t ) and {v,t) on the upper half plane is just the mirror image 
of the path f_(s) between (— fo,^o) an d (— v,t) on the lower half plane. Therefore, in the 
following discussion it is sufficient to consider only a positive initial velocity Vq, without loss 
of generality. Due to this symmetry, the action for the paths v + (s) and f_(s) is the same 
and can be calculated by substituting Eq. ff22l into Eqs. (jUJ) and (jSJ). This yields for the 
basic action 

A[v±,v±]= [ C(v±,v±)ds = A(v,t;v ,t ), (23) 

J to 

where we define 

(e*/^(Ar m + |t;|)-e*o/^(Ar m + t; )) 2 

^ V °> ^ s ' 2r m (e^-e«oA-) " W 

In addition to the set of basic solutions of Eq. (1221) . another formal solution of the EL- 
equation fT2Tl) is given by v(s) = 0. Using the formal solutions v±(s) and v(s) = one 
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FIG. 1: Direct paths in the f-s-plane. (a) A direct path on the upper half plane, parametrized 
by Eq. (|25p . (b) A direct crossing path, parametrized by Eq. (|26p . (c) Here, we plot a number 
of direct paths that are very close together initially. When a direct path crosses the v = axis a 
jump from a positive to a negative curvature occurs, while direct paths that remain on the upper 
half plane continue with a positive curvature. The jump in the curvature leads to a region that 
can not be reached by any direct path from a given initial point (vo,to) (black region). 

can construct the full solution of the EL-equation for fixed initial and final points as linear 
combinations of the three formal solutions. Two distinct classes of solutions then arise, 
namely indirect paths that partly follow the v — axis, and direct paths that do not. They 
are both discussed in the following. 

A. Direct paths 

Direct paths are characterized by continuous v(s) and v(s). They remain either entirely 
on one half of the v-s plane or cross the v = axis. The former are given by the solutions 
v±(s), while paths that cross the v = axis consist of one branch on the upper half plane 
and one on the lower half plane (cf. Fig. Q](a) and (b)). The direct paths (indicated by the 
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subscript d) that remain on the upper plane are simply parametrized by 



\vo,t v +( s )\v ,t ' 



(25) 



where v + (s) is given by Eq. (1221) and the boundary conditions are indicated. 

For the direct crossing paths we have to consider that the upper branch is given by v+(s) 
under the boundary conditions (vo, t ) and (0, t), where i is the time at which the path crosses 
the v = axis, and the lower branch is given by V-(s) under the boundary conditions (0, t) 
and (v,t). Direct crossing paths (indicated by the subscript d and the superscript x) are 
thus parametrized by 



v + (s)\ 



0,t 

vo,to ' 



v -(s)\o} 



t n < s <t 



t <s<t. 



(26) 



Here, the crossing time t is always to < t < t and is determined from the condition of a 
continuous acceleration at the crossover point: 



v + (t) = v.(t). 



Using Eq. (J22J), this condition leads to a fourth order equation for t: 



(27) 



//At„ 



(Ar m + v )e to/Tm + (Ar m - v)e 



t/Tm 



(Ar m + v )e (t0+2t)/Tm + (Ar m - v )e (2to+t)/Tm 



_ Al - (2t +2t)/r m 

where \x = e t//Tm . Eq. ( |28|) has a unique real root t <i < t. 

The action associated with the direct paths that remain on the upper half plane is 

A[v d ,v d ] = A[v+,v+\ = A(v,t;v Q ,t ), 



(28) 



(29) 



which follows immediately from Eqs. (1251) and (|23|) . The action of the direct crossing paths 
on the other hand, consists of contributions from the upper and the lower branch, i.e., 

ri rt 



A[v d ,v, 



/ £(v+, v + )ds + C(v-,v^)ds 

Jto Jt 

A(0,t;v o ,t o )+A(v,t;0,t), 



(30) 



using Eqs. (126|) and (|23|) . Throughout this paper, the function A of Eq. §M§ is adapted to 
the case at hand, by replacing v,t;v , to by the appropriate velocities and times. 
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A crucial observation is then that not all initial and final points in the entire velocity- 
time plane can be connected by a direct path (cf. Fig. 0Jc)). This is due to a jump in the 
curvature of the direct path, when the v — axis is crossed: the upper branch v + (s) is 
always convex with v + rs A close to the v — axis, while the lower branch V-(s) is always 
concave with iL w — A close to the v = axis (cf. Eq. (|2TT) ). There exists thus a region 
in the v-s-plane that cannot be reached by any direct path from a given initial point. As 
t— 1 becomes larger this region grows exponentially, so that eventually, as t — > oo, fewer and 
fewer final points can be reached by a direct optimal path. However, this seems to contradict 
the fact that a stationary distribution exists and has to be approached by f(v, t\vo, to) m the 
asymptotic time limit. Therefore, direct optimal paths can not represent the full solution of 
the EL-equation ( 12TI) . The key is to consider other solutions that satisfy Eq. ( I2TI) piecewise. 
This allows us to construct another class of solutions, namely indirect paths. 



B. Indirect paths 

Indirect paths consist of three parts: a relaxation branch from the initial point (v o, to) to 
the axis at (0,t a ), a part along the zero axis from (0,t a ) to (0, t&), and an excitation branch 
from (0,tb) to the final point (v,t). The relaxation branch is given by v + (s) under the 
boundary conditions (t>o,to) and (0, t a ) and the excitation branch by v + (s) or t>-(s) under 
the boundary conditions (0,fa) and (f,t) (cf. Fig. [2](a) and (b)). Indirect paths (indicated 
by the subscript id) are thus parametrized by 



\0,t a 
I ^0 ,*0 



tn < S < t n 



Vid(S 



{ 



t a < s <t b 



(31) 



Clearly, the times t a , h have to satisfy the conditions t < t a < fa < t. All three parts 
of Vid(s) satisfy the EL-equation piecewise. We note that indirect paths are valid solutions 
of Eq. (|2ip because in the path integral formalism the boundary conditions are the fixed 
initial and final points of the optimal path. If one specifies instead the initial velocity and 
the initial acceleration of the object, indirect paths can not arise. 
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FIG. 2: Indirect paths in the v-s-plane. (a) Two examples of indirect paths on the upper half 
plane, parametrized by Eq. f|3 1 [) . Path 1 (long dashed curve) crosses the v = axis at the times 
t a i and tfji. Path 2 (short dashed curve) at the times t a 2 and t^- (b) Two indirect paths that cross 
the v = axis, parametrized by Eq. (|3ip . (c) and (d): The unique indirect paths with minimal 
action between given initial and final points, parametrized by Eq. (|38p . 



The action associated with these indirect paths is 

rt a rt 

A[v id ,v id ] = / C(v + ,v + )ds + / C(v±, v±)ds 
Jt Jt b 

= A(0,t a ;v ,t )+A(v,t;0,t b ), (32) 

which follows from Eqs. (I3"T|) and (T2"3"|) . Since t a and tb are not specified, there are in principle 
infinitely many indirect paths between (t>o,to) and (v,t) possible, which all have different 
actions according to Eq. (|32l) (cf. Fig. [2]^a) and (b)). But, among all these indirect paths 
there exists a unique indirect path with minimal action. This "optimal" indirect path is 
obtained by determining the minimum of A[vid, v id] of Eq. (|32l) . with respect to t a and tb. 
This is a minimization problem subject to the inequality constraint 



t a < t b - (33) 



In order to find the solution for this minimization we note that 



d d 

— A[v id , v ia ] = ^ A (°> ta, v , t ) (34) 
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is independent of t b and has only one zero for t a G [to? t], at which the axis is crossed with a 
positive slope. Likewise 

d d 

—A[v id ,v id } = —A(v,t;0,t b ), (35) 
ot b ot b 

independent of t a . These properties imply that A[v id , v^] of Eq. ( 1321) has a unique global 
minimum as a function of t a and t b and is monotonically increasing away from it. Due to 
this monotonicity, the minimum of A[v id , v id ], subject to the inequality constraint t a < t b , is 
either given by the global unconstrained minimum or, if this minimum can not be attained 
because it violates the inequality constraint Eq (133]) . by a minimum of A[v id ,v id ] subject to 
the equality constraint t a = t b . 

We then obtain for the minimum of A[vi d , Vi d \. 
(i) The unconstrained minimum is determined by setting each of the time derivatives, 
Eqs. ( 134)) and (|35|) . equal to zero. This yields the times 

t a = t + r m ln (l + -p-) , (36) 
t b = t-T m \n(^l + ^pj . (37) 

We note that the time t a is just the time at which a noise-free (or average) path, described 
by Eq. (jSJ) with £(£) = 0, would relax to the v — axis starting from (v ,t ). Likewise, i b is 
the time at which a noise-free path starting at (v, t) would reach the axis, moving backward 
in time. We parametrize the optimal indirect path specified by i a and i b by (cf. Fig. [2](c) 
and (d)): 

v id (s) = { , i a < s < t b (38) 

, v ±( s )\Zh ' ^ < s < 
The associated action is then given by Eq. (|32|) . where t a and t b are replaced by i a and i b , 
respectively, 

A[d id ,v id ] = A(0, i a ; v , t ) + A(v, t; 0, t b ) 
= A(v,t;0,t b ) 

= U(v). (39) 
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The term A(0, t a ] v , t ) on the right hand side of the first line of Eq. (139]) vanishes, because 
it is the action associated with the noise-free relaxation branch and for these paths the 
Lagrangian is identically zero (cf. Eq. (jHJ)). The third line follows upon substituting Eq. (1371) 
into Eq. I)23j ). 

(ii) If the global minimum does not exist, i.e., if it, < t a , one has to determine the minimum 
of Eq. (|32j) under the equality constraint t a = t b . For this, one has to solve 

d 

— [A(0, t a ] v , t ) + A(v, t- 0, t a )] = 0, (40) 

for t a - One finds that there is a unique real solution of this equation 6 [to , , which is 
identical with the time t determined by Eq. (128]) . The associated action is then given by 
Eq. (132]) with t a = h = t, i.e., A(0, i; vq, t Q ) + A(v, t; 0, t), which is just equal to the action of 
the direct crossing paths, A[v^ ,v^} of Eq. (1311 . 



C. Optimal paths in the v-s-plane 

Having determined the two classes of solutions of the EL-equation (BTI) . one can now find 
the unique optimal path between given initial and final points, direct or indirect, using the 
minimal action principle. To this end we compare the actions of the direct and indirect 
paths separately for the paths on the upper half-plane and for the crossing paths. 

1. We consider the paths on the upper half-plane: The direct paths are given by Eq. (125]) 
with the associated action A[v d ,Vd], Eq. (129]) . The indirect paths are specified by 
Eq. ( |38l) with the associated action A[v i( i,Vid\, Eq. ( |39l . From the condition 

A[v d ,v d ] = A[v id ,v id ], (41) 

which, using Eqs. ( |29l and (139]) . is equivalent to 

U(v) = A(v,t;v ,t ), (42) 

one can derive a critical value of v, denoted by u + (t), such that A[v id ,v id ] < A[v d ,v d ] 
if v < u + (t). Eq. ( ]42l) leads to a quadratic equation for v, which has the relevant root 

u+(t) = e^- t "y T -(Ar m + v )-AT m 

- y/K + 2Ar m v ) (e 2 (t-*o)/r m _ i). (43) 



13 



This means that if v < u + [t) and t > t a , the action of the indirect path is lower than 
that of the direct path. The condition t > i a is necessary for indirect paths to exist. 

2. Let us now consider crossing paths. The direct crossing paths v d (s) are given by 
Eq. fl26|) with the associated action A[v$,Vj], Eq. fl30|) . As before, the indirect paths 
are specified by Eq. (I3"S|) with the associated action A[v id , v id ], Eq. (1591) . The results 
of the minimization of Eq. ( 1321) show the following: When ib > i a the indirect crossing 
path always has a lower action than the direct crossing path between the same initial 
and final points, i.e., A[vid,v i( i} < A[v^,v^]. On the other hand, when t fe < i a , the 
indirect paths of Eq. fl38|) no longer exist and the direct crossing path then has the 
lowest action. Therefore, from the condition 

k = t a (44) 

one can derive a critical value of v on the lower half plane, namely 

u-it) = -Ar m ( A AT ™ e W r <" - l) , (45) 

\AT m +V J 

so that the action of the indirect path is lower than that of the direct path, if v > u~ (t) 
and t > t a . 

It follows from this discussion that for a given initial point {vQ,t ) the optimal path is 
an indirect path if the end point (v, t) lies in the interval u~(t) < v < u + (t) with t > i a , 
otherwise the optimal path is a direct path (cf. Fig. [3]). 

The curves u + (t) and u~(t) represent boundaries in the v-s plane, separating different 
qualitative dynamical behaviors of the moving object in terms of direct and indirect optimal 
paths. In fact, physically, direct paths can be considered as representing a pure slip motion 
of the object, where v ^ apart from one crossing point. Indirect paths follow the v = 
axis for a finite time and thus represent physically a stick-slip motion. 



IV. TRANSITION PROBABILITY 



From the structure of the optimal paths in the f-s-plane, the transition probability 
f(v,t\vo,t ) follows directly via the saddle-point approximation Eq. f|T8|) . The Jacobian 
as given by Eq. (TlOl) reads explicitly for the potential U{v) of Eq. (0J: 

J[ V *] = e (t-to)/(2r m )+Aj* Q S(v*(s))ds^ ^ 
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FIG. 3: Diagram of the optimal paths in a velocity-time plane for a fixed initial point (vo,to) and 
varying final points. The bold black solid curves are given by u + (s), Eq. (|43p . on the upper half 
plane, and by u~(s), Eq. (|45p . on the lower half plane. When the final point of the optimal path 
lies outside of the shaded region, as here (v\,ti), the optimal path is a direct path (solid line). 
Otherwise, if the final point lies inside the shaded region, as here (^2,^2); the optimal path is an 
indirect path (dashed line). 



Here, the term J t * S(v *(s))ds is different from zero only when the optimal path v *(s) crosses 
the v — axis. This means that for the direct paths on the upper half plane the Jacobian 
is just a function of to,t, which can effectively be absorbed into the normalization of the 
transition probability. Only for the direct crossing paths and the indirect paths the Jacobian 
contributes significantly. For the direct crossing paths we obtain 



to 



da = j^jj, (47) 



and for the indirect paths 



6 (v id (s)) ds = 1 + 1 = 1 (48) 
to 2\Vi d (t a )\ 2\v id {t b )\ A 

In deriving these expressions we have used the representation of the delta function [1 

%W) = Er^^-4 (49) 

i \9 \ x i)\ 

where the are the zeros of g(x). In the last step of Eq. (|48j) we have substituted t a and ib 
from Eqs. fl36l) and fl37|) into the time derivative of v id) Eq. fl38l) . 
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With these results for the Jacobian, we can express the transition probability given by 
Eq. ffl8|) as follows. The initial point (v ,t ) is given. 

For t < i a (defined in Eq. (1361) ). no indirect optimal paths occur, so that for a final 
velocity in the interval — oo < v < only the direct crossing paths contribute to the action 
in the expression for the transition probability, Eq. (1181) . while for < v < oo only direct 
paths on the upper plane contribute. Using Eqs. (|29|) and (l30l) . respectively, in Eq. (fl8l) . 
and considering the contribution of the Jacobian, Eq. fj46|) with Eq. fj471) . we obtain the 
transition probability 

-~flA(0,t;v ,t )+A(v,t;0,t)}+A/\v*(t)\ ^ -OO < V < 

(50) 



f(v,t\v ,t ) = A/1 < 



-~fA(v,t;vo,to) 



, < V < OO, 

where M\ is a normalization factor. 

For t > i a indirect paths contribute and the structure of the optimal paths (as discussed 
in Sec. IIII Cj) indicates that direct crossing paths contribute to the action in Eq. (fT8l) for a 
final velocity in the interval — oo < v < u~(t), while indirect paths contribute for u~(t) < 
< u + (t), and direct paths on the upper plane for u + (t) < v < oo. Using Eqs. ( 1291) . ( |30l) 
and (l32i) . respectively, in Eq. ( fl8l) . and considering the contribution of the Jacobian, Eq. ()46l) 
with Eqs. ( H71) and fHHj) . respectively, the transition probability reads 



/(tf,t|u ,*o) =M ( 



-7[A(0,t;u ,to)+A(u,t;0,t)]+A/|^ (t)| 



- 7 f/(t))+l 



-7A(u,t;i;o,to) 



— OO < f < U 



U < f < IT 



M + < f < OO, 



(51) 



where A/2 is a normalization factor. 

Since both u~ — » — 00 and m + — > 00 as t — > 00, only indirect paths contribute to the 
transition probability in the asymptotic time limit and we recover the stationary distribution 
Eq. ([6]) in this limit from Eq. f|5 lj) : 



lim f(v,t\v ,t ) =p(v). 



t—*oa 



(52) 



From the transition probability f(v,t\v ,t ) one can construct joint probability distribu- 
tions for arbitrary sequences of n-points in the velocity-time plane. Due to the Markovian 
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character of Eq. ([I]), the joint probability distribution p(v n ,t n ] ...;v 2 ,t 2 ;vi,t 1 ;vo,t ), which 
contains the probability to find the object at the successive points (v ,t Q ) — > (t>i,ti) — > 
(v 2 ,t 2 ) — > ... — > (v n ,t n ) in the velocity-time plane, is just given by the product of n transi- 
tion probabilities f(v n , t n \v n -i, £ n -i) : 

p(v n ,t n ; ...;v 2 ,t 2 ;v 1 ,ti : v Q ,t ) = f(v n , t„|i> n -i, *n-i) • • • /(v 2 , t 2 \vx, t\)f{v\, ti\v , to), (53) 

where it is assumed that the object is fixed with the initial velocity vq at the initial time to- 
Associated with the joint probability distribution of Eq. (153j) is then an extended optimal 
path along the points (v , t ) — > (v i, ti) — > (v 2 , t 2 ) —>...—»■ (t>„, t n ), which is determined from 
the structure of the optimal paths in the f-s-plane, as discussed in Sec. IIII CI A segment 
of the extended optimal path, between two neighboring points (v n ,t n ) and (f n -i,tn-i), is 
either a direct or an indirect path depending on the relative location of these two points. 

V. CONCLUSION 

We have derived a complete characterization of the optimal paths of the de Gennes' 
equation (CQ) within the path integral framework. The optimal paths can be divided into 
two classes: a) Direct optimal paths, with continuous v(s) and v(s), which physically can be 
considered as representing a pure slip motion of the object, b) Indirect optimal paths, with 
continuous v(s) and discontinuous v(s), which follow partly the v — axis and represent 
physically a stick-slip motion. We have shown that for a given initial point (t>o,to) the 
optimal path will either be a direct or an indirect path depending on the location of the 
final point (v,t) in the velocity-time plane (cf. Fig. [3]). In the asymptotic time limit t — > oo 
a finite final velocity v can only be reached by an indirect path. 

This analysis of the optimal paths leads to an analytical result for the transition prob- 
ability f{v,t\vo,to) in the saddle-point approximation. The calculation of correction terms 
to this result, such as the fluctuation factor F(v*), Eq. ( |20l) . which takes into account the 
second order term in the expansion of the action, Eq. (114H . is left for future work. However, 
we want to emphasize that higher order corrections leave the properties and the structure 
of the optimal paths, and therefore also their slip and stick-slip character, unchanged. 

From a physical point of view, the friction terms in Eq. ([1]) represent a very simple 
phenomenological model of the solid/solid interaction between the object and the surface. 
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By studying generalizations of Eq. (JTJ), incorporating, e.g., two dimensional or memory 
effects, one could model more complicated surface properties, such as surface anisotropies 
or defects. A comparison of such more realistic, but still phenomenological models, with 
experiments will lead to a better understanding of the effects of surface inhomogeneities. 



APPENDIX A: THE PREFACTORS B± AND G± 

The basic solutions v + (s) and V-(s) are given by Eq. (}22l) . where + refers to the upper 
half of the u-s-plane and — to the lower half, respectively. The prefactors B + and C + , for 
the upper half plane, are determined by the boundary conditions, i.e., by the conditions that 
the path is initially at (v ,t ) and ends at (v,t): 

v+(to) = v , v+(t) = v. (Al) 

Solving Eq. fl22|) under these boundary conditions for B + and C + yields 

R e tlTm {v + Ar m ) - e t0 ^(v + Ar m ) 

+ e 2t/T m _ e 2t /T m ' ^ ' 

r e^jvp + Ar m ) - e t0 ^"(v + Ar m ) 

The prefactors B_ and C_, for the lower half plane, are then obtained by simply changing 
A -> -A in Eqs. (El and flA3|) . 
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